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ABSTRACT 

Inhomogeneities in a synchrotron source can severely affect the conclusions 
drawn from observations regarding the source properties. However, their presence 
is not always easy to establish, since several other effects can give rise to similar 
observed characteristics. It is argued that the recently observed broadening of 
the radio spectra and/or light curves in some of the type Ib/c supernovae is a 
direct indication of inhomogeneities. As compared to a homogeneous source, this 
increases the deduced velocity of the forward shock and the observed correlation 
between total energy and shock velocity could in part be due to a varying covering 
factor. The X-ray emission from at least some type Ib/c supernovae is unlikely 
to be synchrotron radiation from an electron distribution accelerated in a non- 
linear shock. Instead it is shown that the observed correlation during the first few 
hundred days between the radio. X-ray and bolometric luminosities indicates that 
the X-ray emission is inverse Compton scattering of the photospheric photons. 
Inhomogeneities are consistent with equipartition between electrons and magnetic 
fields in the optically thin synchrotron emitting regions. 

Subject headings: supernovae: general — stars: mass loss — radiation mecha- 
nisms: non-thermal — scattering 



1. Introduction 



The radio emission observed from supernovae (SNe) can usually be well accounted for 
by synchrotron emission from a homogeneous, spherically symmetric source. In most cases 
though, the limited observations necessitate an assumption about the partition of energy 
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between magnetic fields and relativistic electrons in order to deduce the size and energy con- 
tent of such a source. Although the assumed sphericity finds support in VLBI observations 



of a few nearby radio SNe (e.g., iBrunthaler et al.ll2010l : iBietenholz et al.ll201ll ). the degree 



of homogeneity has been harder to estimate. 

The presence of inhomogeneities would affect the interpretation of observations in at 
least two ways. The low frequency part of the radio spectrum is sometimes seen to deviate 
from that expected from synchrotron self-absorption. This has been ascribed to free- free 
absorption from a thermal population of ele ctrons (e.g., Chevalier l^l982bl: Weiler et al. 2002 ) 
and/or cooling of the relativistic electrons (IFransson &: Bjornssonil998i) . However, some of 
these characteristics could also be caused by an inhomogeneous source structure. Hence, care 
must be taken to include the possibility of inhomogeneities when analyzing such sources. In 
the standard model, the partition of energy between magnetic field and relativistic electrons 
can be determined if either the cooling frequency of the electrons or the inverse Compton 
scattered radiation that they give rise to can be measured. Since the inverse Compton 
scattered radiation is optically thin, it measures the total number of relativistic electrons 
independent of their degree on inhomogeneity. In contrast, the synchrotron self-absorption 
frequency is quite sensitive to the presence of inhomogeneities. Hence, the value deduced 
for the partition of energy between magnetic fields and relativistic electrons depends on the 
degree of inhomogeneity. 



^The winds of hot, massive stars are thought to be driven by line scattering ( ICastor et al. 

I975I). It was soon realized that velocity pe rturbations made them susceptible to an instabil- 
ity (IMacGregor et al.lll979t ICarlberg||l980l ) that would lead to a clumpy wi nd. Several of the 



expected features in such a scenario have been co nfirmed observationally ( iLepine fc Moffat 



I999I : lOskinowa et al.ll2008t ISundqvist et al.ll201ll ). However, although the winds from Wolf- 
Rayet stars, the presumed progenitors of a large fraction of type Ib/c SNe, show signs of 
dumpiness, their main properties do not always conform to those expected from a line driven 
instability ( Oskinova et al.ll2012[ ). The supernova ejecta itself may be clumpy and/or have 
an aspherical dist ribution. The polari zation associated with strong lines in type Ib/c SNe 
has been used by iTanaka et al.l (120121 ) to argue for a clumpy medium. Likewise, double- 
peaked emission li ne profiles in severa l type Ib/c SNe su ggest a flattened distribution of 
the expanding gas (IMazzali et al.ll2005l : iMaeda et al.ll2008[ ). In addition, the shock forming 
process could lead to inhomogeneities; for example, the contact discontinuity is likely to be 
Rayleigh- Taylor unstable. 

The relativistic electrons are thought to be produced by first order Fermi acceleration 
as they repeatedly cross the forward shock front. For a non-relativistic shock, this gives 
rise to an optically thin synchrotron spectral flux F[v) cx with a ~ —1/2 (IBelll Il978 
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Blandford fc Ostrikerlll978l ) or, if synchrotron cooling is important, a ^ —1. Many super- 
novae have a ~ — 1. For a supernova exploding into the wind of a red super giant, the 
energy d ensity behind the shock caii be high enough for synchrotron cooling to become im- 
portant (IFransson &: Bjornssonlll998l ). However, for Wolf-Rayet stars the density of the wind 
is thought to be so low that such cooling is unlikely. Inverse C ompton cooling may be signif- 
icant but only during the very early phases (e.g., SN 2002ap; iBjornsson fc FranssonI 12004 ) . 
A possible way to reconcile the observed spectra in type Ib/c SNe with a diffusive shock 
acceleration origin for the relativis tic electrons is tha t the acceleration is so efficient that the 



whole process becomes non-linear (IBaring et al.lll999l : iBerezhko fc Ellisonlll999l : lEUison et al. 



20001 ) . This results in a concave spectrum. The steepening at low frequencies, which is rele- 
vant for the radio regime, is determined by a combination of the injection efficiency and the 
escape of the highest energy electrons from the acceleration process. Neither of these are 
well understood physically. Hence, if it could be established that the acceleration process 
in type Ib/c SNe is non- linear, the observed value of a ~ — 1 would strongly constrain the 
physical mechanisms operating in the shock region. 

Inhomogeneities may provide an alternative to a non-linear acceleration process as the 
origin for the observed steep spectral flux. In the standard scenario they could modify the 
acceleration process so that a steeper distribution of electron energies results. On the other 
hand, if the acceleration process is associated with the contact discontinuity/the reverse 
shock region, the inhomogeneities may also here leave some imprint on the observed spectral 
flux. 

A rel ated issue concerns the origin of the X-ray emission in SNe. This is still a matter 
of debate (IDwarkadas fc Gruszkdl2012l ): for example, in SNe with detectable radio emission 
a straightforward extrapol ation of the synchrotron spe ctrum falls substantially below the 
observed X-ray luminosity. IChevalier fc FranssonI (120061 ) argued against a thermal origin in 
type Ib/c SNe. Instead they showed that with a suitable choice of parameters both the 
radio and the X-ray emission can be attributed to synchrotron radiation from a concave 
distribution of electron energies expected from a non-linear acceleration process. Another 
alternative is inverse Compton scattering of the photospheric photons. Such an origin for 
the X-ray emission in SN 2002ap is co nsistent with a homogeneous, spherically symmet- 
ric source ( Bjornsson &: FranssonI 2004). However, SN 2002ap stands out in that its radio 
flux peaked early (within a few weeks). For SNe peaking at later times, the radius of the 
shock tends to be larger and the optical emission has started to decline. In such SNe, 
homogeneous and sphericall y symmetric models usually under predict the X-ray emission 
(IChevalier fc FranssonI l2006l ). This conclusion hinges on the assumption that the energy is 
partitioned roughly equally between relativistic electrons and magnetic flelds. 
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As mentioned above, the partition of energy and the degree of inhomogeneity are in- 
tertwined. The aim of the present paper is to discuss how the observed radio emission can 
be used to infer the presence of inhomogeneities and the imphcations it has for the X-ray 
emission. Focus is on acceleration at the forward shock. However, the conclusions are general 
enough to be applicable with only minor changes also to the case when the emission comes 
from the contact discontinuity/the reverse shock region. The scope of the paper is limited 
to type Ib/c SNe. The properties of an inhomogeneous synchrotron source are discussed in 
§2. Special attention is given to characteristics that are unlikely to be produced by other 
mechanisms. The non-linear diffusive shock acceleration scenario in a supernova context is 
discussed in § 3. It is shown that it has a few implications, which may constrain its ap- 
plicability. The consequences of an inhomogeneous source structure for the observed X-ray 
emission are considered in § 4, assuming the latter to be produced via inverse Compton scat- 
tering of the photospheric photons. It is shown how the combined effects of the partition of 
energy and inhomogeneities can be deduced from the observed radio, bolometric and X-ray 
luminosities. A discussion of the results as well as a brief summary of the main conclusions 
follow in § 5. It is pointed out that the presence of inhomogeneities can have implications for 
the relation between ordinary Ib/c SNe and engine-driven hypernovae as well as the basic 
properties of the acceleration of relativistic electrons. 



2. An inhomogeneous synchrotron source 



With the increased quantity and quality of radio observations of SNe, it has become 
clear that spectra and/or light curves are often broader than expec ted for a standard ho- 
mogen eous, spherically symmetric source. In the type lib SN 1993J, iFransson fc Bjornsson 
(Il998l ) modeled the broad light curves as due to cooling of the synchrotron emitting elec- 
trons. Such an explanation requires a high density behind the shock. Since the wind velocity 
of the progenitor stars to type Ib/c SNe are expected to be much higher than those for the 
progenitor stars to type lib SNe, the densities will be correspondingly smaller. Hence, the 
broad spectra/light curves in type I b/c SNe are not likely to be caused by cooling of the 
same magnitude as that deduced by IFransson fc BjornssonI (119981 ) in SN 1993J. 



Soderberg et al.l (120051) introduced a parameter to artificially broaden the spectra in the 



type Ib/c SN 2003L. IWeiler et al.l (1201 ll ) discussed various causes for the spectral shape in 
the type Ic SN 19941. They concluded that the best fit to the observations came from free-free 
absorption intrinsic to the synchrotron source. The thermal electrons need to be distributed 
roughly as the relativistic ones, since, in order to get a good fit, the radio opacity should 
be the sum of synchrotron and free-free absorption. Since there are no clear indications of 
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extrinsic free-free absorption in type Ib/c SNe, the thermal gas needs to cool after passing 
through the shock. With the circumstellar densities expected in type Ib/c SNe, this suggests 
the presence of dense clumps of t hermal gas. For a clumpy gas the need for cooling sets 



a lower limit to their typical size f lChugailll993l ). The effects on the synchrotron spectrum 
will be quite different depending on whether such clumps are optically thick or not to free- 
free absorption. Optically thick clumps produce a depression of the synchrotron spectrum 
without substantially c hanging its shape, except at the higher frequencies where the clumps 
become optically thin (jChugai fc Chevalierll2006l ). Hence, if the observed broadening of the 
radio spectra in SN 19941 is due to internal f ree-free a bsorp tion, it implies optically thin 
clumps. However, following the arguments in IChugail ( 119931 ). the combined requirements 
of cooling and optically thin clumps are unlikely to be met in standard type Ib/c SNe 
environments. 

Below, deviations from a homogeneous, spherically symmetric source will be described 
in terms of the covering factor and the filling factor. The former characterizes the optically 
thick properties of the source while the latter accounts for its optically thin properties. One 
important aspect that does not enter in an analysis of an unresolved source is the geometrical 
distribution on the inhomogeneities. Hence, in this paper, the concept of inhomogeneities 
includes also an aspherical distribution of emitting plasma. 



2.1. Inhomogeneities with a covering factor less than unity 

The spectral distribution of the optically thick emission in an inhomogeneous source 
is the combination of two independent factors. Since the synchrotron source function is 
proportional to B^^^'^, where B is the magnetic field strength, the variation of the covering 
factor /b,cov with B needs to be specified. Furthermore, Fl^u) is proportional to the brightness 
temperature. At a given frequency, this, in turn, is given by the value of the Lorentz factor 
(7abs) corresponding to a synchrotron optical depth of unity. Since 7abs c<: {u/By^"^, the 
second relation that needs to be specified is that between B and u. 

The source is assumed to consist of a spherical shell of outer radius R and thickness R/A. 
Let the observed emission along a given line of sight through the source be dominated by a 
magnetic field B over a distance r{B) < R/A. The probability to find a particular value of B 
is taken to be P{B) = A B^"^ for Bi < B < B2 and zero outside this range. Here, A oc B^^^ 
for a > 1. The spectrum has three different regimes: For u < u{Bi), the spectrum is that 
for an optically thick source with magnetic field B = Bi, while for u > u{B2), the spectrum 
corresponds to an optically thin source with magnetic field B = B2. In the transition region 
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z/(i?i) < u < i'{B2), the spectrum is given by 



F{iy) oc R^iy 



B2 

2„5/2 / PiB) 



dB. 



(1) 



B 



Hence, F{u) oc R^u^^^B^-^B^/^-". 



The relation depends on the more detailed properties of the source. The distri- 

bution of Lorentz factors (7) for the relat ivistic electrons is taken to be A^(7) oc 7"^ for 
7 > 7min and p > 2. Since observations (IChevalier fc FranssonI |2006| ) indicate p ~ 3, in 
order to simplify the discussion, p = 3 is used in the following. In the transition region, 
z/^ oc ?7e%'"(-S)7min/7abs- Here, Ue and Ub = -B^/Svr are, respectively, the energy densities 
in relativistic electrons and magnetic fields. Together with the expression for 7abs, this gives 



B oc 



7/5 



(f/er(5)7^i,)'/^' 

The simplest case is when only the value of B varies so that B oc u"^^^, which leads to 



(2) 



oc z/ 



(16-7a)/5 



(3) 



The observations of SN 2003L by lSoderberg et al.l (120051 ) indicate that the broadening of the 
spectral peak is quite symmetric implying a rather flat spectrum in the transition region. 
This is obtained for a ~ 16/7. The corresponding covering factor is /b,cov ~ (B/BiY^"-, 
which, together with equation ([2]), gives a covering factor varying with frequency as fu,cov ~ 
iu/uiB^))-^/\ 



The parametrization of the spectra for SN 19941 done by IWeiler et al.l ( 120111 ). can be 
used to estimate the frequency width of the transition region. The parameter values derived 
by them indicate that this width slowly decreases with time but is typically around two in 
the radio regime; i.e., v{B2)/i'{Bi) ^ 2. Hence, the covering factor for the optically thin 
synchrotron emission (i.e., v > i^{B2)) is roughly 1/4. Although ISoderberg et al.l ( l2005l ) used 
a different parameterization to account for the broadening of the spectral peak, the implied 
value for the covering factor is consistent with that deduced for SN 19941. 



2.1.1. Implications from the light curves 



Light curves provide complementary information to spectral data. However, an extra 
assumption has to be made regarding the mass loss rate and wind velocity of the progenitor 
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star. The discussion in this section will assume that both of these are constant. Although this 
may be a fair approximat ion for the majority of type Ib/c SNe, t here are possible exceptions ; 
for example, SN 2003bg JSoderberg et alibooeh and SN 2001em Jchugai fc Chevaheiibooel ). 



The three different spectral regimes translate into three different regimes for the light 
curve. It is seen from equation ([1]) that in the transition regime, the flux varies with time t 



as F(t) cx R B1 . It is usually assumed that the energy density of the magnetic field scales 



with t he thermal energy density behind the shock, which implies B oct ^. iMatzner fc McKee 



( I1999I ) showed that the high shock velocities rele vant for type Ib/ c SNe is expected to scale 



with time roughly as ^ hence, F{t) cx t'^-^'"'. ISoderberg et al.l (120051 ) find f or SN 2003L 
that (X t^'^ before the peak. A similar behavior can be seen in SN 19941 ( iWeiler et al. 



2OIII ). This results in a value of a somewhat smaller than deduced in section [2?T] from the 
spectral shape. It is seen from equation ([T]) that this implies a more rapid variation of the 
magnetic field strength with frequency than obtained with the above assumptions. 

The relation is determined by equation ([2]), which shows that a steeper variation 
will result from a decrease of r{B) and/or Ue with increasing value of B. Since the covering 
factor decreases with the strength of the magnetic field (a > 1), it is reasonable to expect 
that the same should apply for r{B). On the other hand, a decreasing value of Ue with B 
is non-trivial and constrains the mechanism(s) that accelerates the electrons and amplifies 
the magnetic field. An interesting case is when the total energy density/pressure of the 
relativistic electrons and magnetic field is constant behind the shock. Since, in the transition 
zone, u oc (f/gt/g)^/^, the conditions in the regions producing the majority of the optically thin 
radiation (i.e., B = B2) correspond, roughly, to equipartition between relativistic electrons 
and magnetic fields. The implications of this, with respect to the origin of the observed 
X-ray radiation, will be discussed further in section HI 

A smaller value of a will affect the covering factor (/b.cov ~ {B/BiY^""). However, as 
discussed above, the implied range of 5- values in the transition regime increases so that the 
value of /b2,cov is only marginally smaller. 



2.2. Inhomogeneities with a covering factor of unity 

The type of inhomogeneities discussed in section 12.11 is consistent with a scenario in 
which the wind of the progenitor and/or the ejecta of the supernova is clumpy. There is 
another type of inhomogeneity where the source is spherically symmetric (i.e., the covering 
factor is unity) but the values of the energy densities of the magnetic field and relativistic 
electrons both vary systematically with distance (r) away from the shock front. In this 
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case, the inhomogeneity is likely due to the processes that amplify the magnetic field and 
accelerate the electrons. The optically thick spectral index is determined by the different 
variations behind the shock front of the synchrotron emissivity and absorptivity. Since their 
relative variation (i.e., the source function) depends on B only, the main features can be 
illustrated by assuming a constant density of electrons. Such a situation can be described 
by 



B 



Bo: 
Bo 



r < Tn 



To <r < R/A 



(4) 



Let be the value of r corresponding to an optical depth of unity at the frequency u. 
Then for v < z/q, where r^^ = Tq, the spectrum is that for an optically thick, homogeneous 
source. Likewise, for frequencies such that > R/A, corresponding to z/ > i/r, the syn- 
chrotron emission is that for an optically thin source. Hence, in the same way as in section 
12. 1[ there is a transition region between the standard optically thick spectrum corresponding 
to a homogeneous source and the optically thin radiation. 



The synchrotron optical depth r is 



r oc 



r ^5/2 
7/2 



dr. 



(5) 



With the use of equation (jl]), in the frequency range z/q < z/ < z/r, this leads to 



z/ oc < 1 + 



5b/2 + 1 



5b/2+l 



2/7 



(6) 



Likewise, the spectral flux in the transition region is given by 

'0 



oc 



-dr 



1 



oc u 



2b + 1 



2b+l 



(7) 



In the limit when the contributions to both the emission and absorption increases with 
distance away from the shock front (i.e., b > —2/5), equations ([6]) and (171) can be combined 
to give 



a 



5 + 9b 



2 + 5b 

It is seen that for b > 0, the spectrum is flatter than for the homogeneous case. However, the 
smallest value is a = 9/5, which, as discussed above, is too large to account for the observed 
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spectra. In the other hmit (i.e., b < —1/2) the main contributions to both emission and 
absorption come from r ^ r^. Although this can give rise to flat spectra (a ~ 0), it can be 
deduced from equations dHD and ((Tj) that the spectral range of the transition region is too 
small to accomodate the observations. Since, in the range — l/2<b< — 2/5, the spectra 
rise more steeply (a > 5/2) than for the homogeneous case, it is unlikely that the type of 
inhomogeneities discussed here are the cause for the broadening of the observed spectra and 
light curves in type Ib/c SNe. 



3. X-ray emission due to synchrotron radiation from a modified electron 

distribution 

The repeated scatterings across a shock, inherent to diffusive shock acceleration of parti- 
cles, impart momentum to the thermal gas in the upstream region. In order to conserve mo- 
mentum and energy, this pre-acceleration of the upstream gas implies tha t energy/particles 



must escape from the flow/acceleration process (IBerezhko fc Ellisonlll999l ). When the injec- 
tion efficiency is low, the effects on the upstream gas are negligible and the flow is character- 
ized only by the shock compression ratio. At some critical value of the injection efficiency, 
the upstream flow is signiflcantly affected by the acceleration process and particle escape 
needs to be taken into account. In addition to the shock compression ratio (rgub)? the flow 
is now characterized also by the total compression ratio (rtot)- In this case, the flow is de- 
termined not only by the injection efficiency but also the momentum above which particles 
escape the acceleration process. Likewise, the energy distribution o f the accelerated par ticles 



is modified by the changing flow characteristics. As discussed by Ellison et al.l ( l2000l ). the 
standard energy distribution (i.e., p = 2) is roughly valid in an intermediate energy range 
but steepens at low energies {piow > 2) and flattens at high energies (phigh < 2). 

The value of pugh is rather model independent as long as the acceleration is in the non- 
linear regime. The flow upstream of the shock is then determined by the momentum input 
from the high energy part icles. A self consistent description of this process gives Phigh ~ 1.5 



(IBerezhko fc Ellisorull999l ). The value of piow is related to rsub in the same way as in standard 
shock acceleration. However, in non-linear shock acceleration the value of rsub results from 
a self consistent calculation of the escaping energy/particles. Hence, the value of piow is 
more model dependent than that of Phigh, since it is sensitive to the details of both the 
injection process and the escape of particles. Neither of these processes are well understood 
physically and are normally treated by introducing two or more parameters. Furthermore, 
the calculated energy distribution is that for the momentum carrying particles (i.e., ions). It 
is usually assumed that the diffusion length depends only on energy (e.g., Bohm diffusion). 
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so that the energy distribution of electrons can be taken to be similar to that of the ions. 
The type Ib/c SNe studied by I Chevalier &: FranssonI (120061 ) have piow ~ 3 (i.e., a ~ — 1) 
with a rather small spread, which corresponds to r^uh ~ 2.5. In contrast to the generic value 
Phigh ~ 1-5, this value of piow constrains the physics of non-linear shock acceleration. 

The energy density of the magnetic field as well as that of the low energy part of 
the electron distribution (i.e., the energy density at the injection momentum) are normally 
assumed to be some (constant) fraction of the thermal energy density behind the shock. 
Although the total compression ratio is larger for non-linear as compared to standard shocks 
(i.e., rtot > 4), the thermal energy density behind the shock is smaller. Hence, the implied 
mass loss rate of the progenitor star increases correspondingly. For rsub ~ 2.5, standard shock 
relations show that the thermal energy density behind the shock decreases by a factor rtot- 
Furthermore, the influx of energy decreases by a factor rtot/4, decreasing the column density 
of magnetic field and relativistic electrons by the same amount. This increases the implied 
energy density of the magnetic field by a factor (rtot/4)^/^^ {p\ow = 3). Taken together, these 
effects increase the needed mass los s rate of the progenitor star by a factor 4(rtot/4)^^''^^. 
The numerical calculations done by I Chevalier &: FranssonI (120061 ) to fit the observations of 
SN 19941 have '"tot ~ 46. The implied mass loss rate of the progenitor star is then more 
than a factor 10^ larger for non-linear shock acceleration as compared to standard shock 
acceleration. 

It could be argued that the energy density of the magnetic field scales with the total 
energy density behind the shock rather than the thermal one. In this case Ub ~ '"totf^e"™? 
where Ul°'^ is the energy density of the relativistic electrons corresponding to the piow-part 
of the energy distribution. The implied value for the magnetic energy density would increase 



by a factor 
then : 



„8/19 
' tot 



tot 



(Plow = 3). The needed increase of the mass loss rate of the progenitor is 
/4)*/^^, which is roughly a factor 10 lower than for the case in which the energy 



density of the magnetic field is assumed to scale with the thermal energy density. 

Adiabatic heating may not be the only process which determines the temperature struc- 
ture upstream of the shock. Additional heating pro cesses will decrease the co mpression ahead 
of the shock and, hence, decrease the value of r tnt (IBerezhko fc Ellisonlll999l ). It is seen from 
the calculations by I Chevalier &: FranssonI (120061 ) that Phigh ~ 1-5 is needed to simultaneous fit 
the radio and X-ray emission from SN 19941; basically, a flattening at higher energies of the 
electron distribution is required to offset the steepening of the spectrum due to synchrotron 
cooling. In order to keep the generic Phigh ~ 1-5, the value of net cannot decrease too much . 
The estimates of the mass loss rates of the progenitor stars in I Chevalier fc FranssonI (120061 ) 
was done using radio data interpreted within a standard shock scenario. Since their deduced 
values are already quite high, assuming the presence of non-linear shock acceleration may 
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well require prohibitively large mass loss rates. 



Since the electrons radiating in the X-ray regime are cooling, the time evo lution of the X- 

ray lu minosity {L^) from a homogeneous, spherically symmetric source is (e.g.. lBjornsson &: Fransson 
200l) 

oc viB^^-'y\ (9) 

where Vsh is the shock velocity. Assuming a constant mass loss rate and that the energy 
density of the magnetic field scales with the thermal energy density b ehind the shock, B oc 



t ^. Furthermore, in self-similar solutions for the shock propagat ion (IChevalierlll982al). th e 



high velocities relevant for type Ib/c SNe lead to Vsh oc t~°-^^ ( Matzner fc McKe3 1999 ). 
With p ^ Pugh ~ 1-5, one finds L^'^t'^'^^. Hence, the X- ray luminosity is expe c ted t o 
decline slowly. This is consistent with the numerical results in I Chevalier fc Fransson! ( l2006l ): 
for example, in their calculations the X-ray luminosity of SN 19941 declines by less than a 
factor 2 from t = 100 days to t = 1000 days. Except for an apparent increase in SN 19941 
during its initial phase, the observed X-ray variability in SNe 19941, 200 lig, and 2003bg 
( jChevalier &: FranssonI 120061 ) suggests a decrease considerably faster than predicted from 
equation (El). This is i n line with the observed X-ray variability in other typ e Ib/c SNe (e.g., 
Soderberg et al.l 120081 ). It was pointed out by I Chevalier fc FranssonI ( 120061 ) that the X-ray 
emission in SN 19941 on days 82 and ~ 2500 could not both be due to synchrotron radiation 
from electrons accelerated in a non-linear shock. 

Taken together, the issues raised above show that invoking non-linear shock acceleration 
to explain the X-ray emission meets with some complications. Therefore, in the next section, 
the merits of inverse Compton scattering of the photospheric photons as the origin of the 
X-ray emission are considered. 



4. X-ray radiation from inverse Compton scattering of the photospheric 

photons 

During the first few weeks after explosion, inverse Compton scattering of photospheric 
photons can contribute significantly to the X-ray emission in type Ib/c SNe. At later times, 
such an origin for the observed X-ray emission is more model dependent; for example, a 
homogeneous, spherically symmetric source, in which the energy densities in magnetic fields 
and relativistic electrons are in equipartition, falls short by at least an order of magnitude 



( IChevalier fc FranssonI l2006l ). However, in light of the discussion in section |2l an inverse 
Compton origin of the X-ray radiation may still be viable. In this section, the assumption 
of homogeneity is relaxed with the aim to elucidate how this affects the deduced source 
properties. The type of inhomogeneities considered is limited to those discussed in section 
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EH 

Let Lj. and Lboi denote the optically thin radio luminosity and the bolometric luminosity 
of the supernova, respectively. Since the optically thin radio spectral index is close to a = — 1, 
Lr cx vFiy) is roughly independent of frequency. In order to facilitate the discussion it will 
be assumed that the distribution of relativistic electrons has the same form throughout the 
source although its amplitude may vary. 

Since the inverse Compton scattered radiation is optically thin, it is unaffected by 
inhomogeneities and depends only on the total number of relativistic electrons. With < Ue > 
denoting the average energy density of relativistic electrons, the inverse Compton scattered 
radiation can be written Lx/Lboi = (4/3)7^ r(7), where r(7) is the scattering optical depth 
of electrons with Lorentz factor 7. Hence, 

2 , , <?7e>7minP-2 / 7 V"^ R 

7 ^(7) = 2 7 (^T-r 

= ^^^^^^Ti? forp = 3, (10) 

where ctt is the Thomson cross-section. It has been assumed that the thickness of the shell 
is i?/4 and p = 3 corresponds to a = — 1. The average energy density of electrons can then 
be expressed as 

6mc^ Lx , , 

< >= (11) 

7min<7T-K -^^bol 

The synchrotron emitting regions that contribute most of the optically thin radio emis- 
sion is assumed to correspond to a volume filling factor /b,vo1- In order to simplify the 
notation, in this section B corresponds to B2 in section 12.11 so that, for example, /b,vo1 ~ 
/b2,cov'"(-B2)/(-R/4). The energy density of the magnetic field in these regions is denoted 
by Ub- Since synchrotron radiation can be regarded as inverse Compton scattering of the 
virtual photons due to the magnetic field, 

Lr _ /b,vo1 Ub^ ^^2') 

Here, f/ph = Lboi/47ri?^c is the energy density of photons due to the photospheric emission 
and S = Ue/ < Ue >, where Ue is the energy density of electrons in the synchrotron emitting 
regions. The partition of the energy densities of relativistic electrons and magnetic fields in 
the regions responsible for the optically thin radio emission is then given by 

-/b,voi5 - — 1- 7 — (13) 



Ub O'T ' 7min-f^r V-^ 



^bol 
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4.1. Mass-loss rates and shock velocities 

In order to minimize the mass-loss rate of the progenitor star it is often assumed that all 
the available electrons are injected into the acceleration process. This allows a determination 
of both the mass-loss rate and the value of 7min- With the assumption that a fraction ee of 
the thermal energy density behind the shock is converted into relativistic electrons, 

9 p-2 rup /Vshy 



7mm = ^ rCc /ie, (14) 

where rup and He denote, respectively, the proton mass and the mean molecular weight of 
the electron and, for p = 3, this can be written 

7mm = 1.9 X 10(3 ee)^f^V. (15) 



2 \ c ^ 

For a progenitor wind devoid of hydrogen and complete ionization behind the shock /Xe ~ 2 
The mass-loss rate (M) is related to < f/c > through 



9eeM/^;3h\2 



where is the wind velocity of the progenitor. With the use of equation ffTTj) . this implies 



l.SxlOM— -^T^. (17) 

fw,3 \3v^hJ (3ee)^bol 

Here M_5 = 1 x 10~^ Mq yr~^, Ww,3 = fw/(10^ km s~^), and the observing time t is written 
as ti = t/lOdays. Furthermore, the expression for 7jnin in equation f[T^ has been used. 

In an inhomogeneous source, the shock radius is not directly related to the radiating 
volume. Since only a fraction (i.e.,/B,voi) of the volume behind the shock front contributes, 
the shock velocity needs to be larger than for the homogeneous case. Again, in order to 
simplify the notation, = ^^(-82) will be used. Standard synchrotron theory then gives for 
a spherical source. 



'^mJB,cov / 



FiVr^yl^' Vf^e(7abs)A^ 



B oc 



where ?7e(7abs) is the energy density of relativistic electrons radiating at frequencies larger 
than i/m and A = Ar{B)/R. Using the expression for < f/c > in equation ( fTTI) . equations 
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f fTSjl and f lT9|) can be written 



B « ^^-(^^)'^ (20) 
R « Q/in „i/2 • (21) 



9/10 .1/2 \ A 5 

Since Lr ~ ^'m-^('^m), it is seen that equation ( |T3l) is implied by equations ( 120|) and lpTl) . 



4.2. Comparing SN 2003L to SN 2002ap 



Although several of the unknown parameters characterizing the inhomogeneities affect 
the value of the radius and, hence, the shock velocity, it is seen from equation fl2Tl) that the 
main effect is due to the covering factor so that fshSS/Bcov- discussed in section [2?T| 
the value of /b,cov can be estimated from the broadening of the radio spectra and/or light 
curves. The radio luminosities of type Ib/c SNe span a rather wide range with SN 2002ap 
(low) and SN 2003L (high) as extreme cases. Although the observations of SN 2002ap are 
not as extensive as those of SN 2003L, there are no indications of spectral broadening. Since 
/b,cov ~ 1/4 is deduced for SN 2003L in section [2711 its actual shock velocity would then be 
a factor ~ 2 larger than deduced for the homogeneou s case. This makes the velocity almost 
equal to that derived in I Chevalier fc Franssoru (120061 ) for SN 2002ap. 



The volume filling factor of the optically thin synchrotron emission can be estimated 
from equation f[T^ . 



Ub 



-^ = 2.8x10Vb,voi5'^'° 



360 317811 \-^^bol,L/ -^^r,L 



(22) 



The observed values are normalized to those given by ISoderberg et al.l (120051 ) for SN 2003L 
at t 40 days (i.e., = 9.2 x 10^^ erg s-\L, = 0.074 x = 6.8 x 10^^ erg s-^ and 
-^boi = 1.5 X lO''^ erg s^^). The value of /b,vo1 depends on the relative distribution of magnetic 
fields and relativistic electrons and, hence, it is harder to estimate than that of /b,cov from 
the available observations. However, as was discussed in section 12. H the combination of light 
curves and spectra can be used to constrain the source properties. 

One way to account for the radio light curves in SN 2003L is to assume an anti- 
correlation between the energy densities of electrons and magnetic fields. In order to il- 
lustrate how this can be used to deduce the value of /b.voI, consider the special case when 
their sum is constant throughout the source. This implies that the optically thin synchrotron 
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radiation comes mainly from regions characterized by rough equipartition between electrons 
and magnetic fields. Since 5 1/2 for equipartition, one finds from equation fl22|) that 
/b,vo1 ~ 1/70 for Cc = 1/3. Furthermore, with /b,cov ~ 1/4, the thickness of the emission 
region is A ^ 1/17. 



Comparing SN 2002ap to SN 2003L, one finds from the values given in I Chevalier fc Fransson 



( 120061 ) that at t ^ 4 days (i.e., a factor 10 earlier than for SN 2003L), (Lx/-Lboi)^/-^^r is a 
factor two smaller. Hence, the RHS of equation fl22p is a factor 20 smaller for SN 2002ap 
as compared to SN 2003L. Assuming, for example, equipartition between relativistic elec- 
trons and magn etic fields in the optica l ly thi n emission regions then implies /b.voI^^ ~ 0.1. 



The reason that iBjornsson fc FranssonI (j2004[ ) found that the observations of SN 2002ap are 



consistent with a homogeneous source in equipartition is that they assumed cooling to be 
important. Since, in this case p ~ 2, the energy density in relativistic electrons varies only 
logarithmically with 7min- On the other hand, for p = 3, the value of 7min is important since 
then Uc oc l/7min- Hence, assuming cooling to be negligible requires either Uc/Ub ~ 10 or 
an inhomogeneous emission region also in SN 2002ap. However, in the latter case, the inho- 
mogeneities are likely to be less pronounced than in SN 2003L; in particular, observations 
suggest /b,cov ~ 1 in SN 2002ap. 



4-4 U:- TT^- (23) 



The mass-loss rate derived from equation f lT7|) for the progenitor star to SN 2003L on 
t = 40 days is 

Vw,3 y^VshJ (3ee)' 

This is roughly a factor 10^ larger than the corresponding mass-loss rate for SN 2002ap. 
Hence, interpreting the X-ray emission as inverse Compton scattering of photospheric pho- 
tons leads to the conclusion that the main difference between SN 2003L and SN 2002ap is 
the much larger mass-loss rate associated with the former supernova. 



4.3. Implications for other SNe of type Ib/c 

Although SN 2003L and SN 2002ap have rather different observed properties, the value 
of (Lx/Lboi)^/-f'r differ only by a factor two. Figured] shows that variations in this quantity 
are quite small also amongst other type Ib/c SNe; in particular, in comparison with the large 
variations of the radio luminosity. The exception is the late time observations of SN 19941. 
The observations at t ~ 2500 days give a very large value and are not included in Figured] 
This tendency of an increasing value of {L^/ Ly^oiY / can be seen already at t ~ 82 days. A 
possible explanation is an additional source of X-ray emission, which is supported by the fact 
that the X-ray luminosity of SN 19941 actually increased between days t ^ 52 and t ^ 82. 
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It should also be noticed that a certain amount of scatter in (Lx/-Lboi)^/-Z^r is expected since: 
(1) Different energy bands were used to calculate Lx- (2) In the derivation of equation 
V = 3 was assurned. T he radio spectra indicate some variation around this value and 



Chevalier fc FranssonI (120061 ) used the observed spectral index to deduce at 5 GHz. 



Rewriting equation (!22l) or using equations ( 120|) and (pT]) leads to 



^bol 



^ 1 UeeeVsh 1 / ,x 



Although the observations are still rather sparse, the small scatter in Figured] shows that 
they are consistent with L^%{L^/ L^o\Y during the first few hundred days of interaction 
between the supernova shockwave and the circumstellar medium. This is as expected for 
a scenario in which a substantial contribution to the X-ray emission comes from inverse 
Compton scattering of photospheric photons. 

For a homogeneous source the partition of energy between relativistic electrons and 
magnetic fields is directly obtained from {L^/L]^oiY / L^. Although its value is not so straight- 
forward to deduce for an inhomogeneous source, equation ( l24l) indicates that one may expect 
the value of (Lx/-Lboi)^/-^r to decrease with time for an individual source. Unfortunately, 
X-ray observations of individual sources are not yet extended enough for such a comparison 
to be made. However, it is seen in Figure[2]that for the SNe discussed in this paper, the value 
of (Lx///boi)^/-^r vary only weakly, if at all, with time. Unless these SNe happened to be 
observed at particular times, this suggests that the same conclusion holds also for individual 
sources. Hence, an inverse Compton scenario for the X-ray emission implies not only an 
inhomogeneous source structure but also that variation of the inhomogeneities in type Ib/c 
SNe to some extent cancels variation in t. 

To some extent the properties of the inhomogeneities can be constrained from the de- 
tailed radio observations. The optically thin radio luminosity Lj. cx 7minf4f^B-R^/B,voi- With 
the st andard assumption that U p oc Ub oc 1/t^, this leads to Lj. oc v L fB,m\/t. With oc 



(IMatzner fc McKed Il999l ). the observed time variation of L,- (j Chevalier fc Fransson 
20061 ) indicates /b,vo1 ~ constant (SN 2002ap is an exception; its radio luminosity decreases 
considerably slower than those for the rest of the type Ib/c SNe). With a roughly constant 
value of /b,vo1, equation (IMl) suggests that the main variation with time is a decreasing value 
of 6, i.e., the average energy density of relativistic electrons (< Ue >) decreases slower with 
time than that in the synchrotron emitting regions {Ue). 
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5. Discussion and conclusions 

Compared to a homogeneous synchrotron source, the presence of inhomogeneities im- 
phes a higher velocity for the forward shock. The observed radio spectral index [a ~ —1) 
makes the injection energy of the relativistic electrons/low energy cut-off (7min) an impor- 
tant parameter, since the total number of relativistic electrons oc 7~fn- Together with the 
assumption that all the electrons are injected into the acceleration process, inhomogeneities 
lead to a lower mass-loss rate for the progenitor star. On the other hand, the higher velocity 
of the forward shock increases the total energy (and/or decreases the mass of the ejecta). 
Hence, if the physical parameters of a source are deduced using a homogeneous model, a 
varying covering factor (/b,cov) would result in an apparent correlation between total energy 
and shock velocity. This could, to some extent, account for the difference between ordinary 
type Ib/c SNe and those suggested to be powered by a central engine. 

The value of /b.cov deduced from the broadening of the spectra and/or light curves gives 
a maximum value, since it is possible that not the whole source contributes to the optically 
thick synchrotron radiation. With this maximum value of /b,cov, it was shown that the 
velocity of the shock in SN 2003L would be similar to that in SN 2002ap. However, smaller 
values of /b,cov and, hence, higher shock velocities cannot be excluded. The SNe in Figured] 
are all ordinary type I b/c except for SN 2009bb, which is suggested to have an engine by 



Soderberg et al.l (120101 ). It is seen that SN 2009bb follows the same relation as that of the 
ordinary SNe; in particular, its radio luminosity does not stand out in comparison with 
those of the most luminous of the ordinary type Ib/c SNe. The main difference is instead its 
low value of the synchrotron self-absorption frequency {u^) at comparable epochs. For an 
electron energy distribution with p = 3, the optically thin radio luminosity oc t'j^/B.cov 
It is therefore possible that the low value of the synchrotron self-absorption frequency is 
compensated for by a larger covering factor. This then necessitates an even higher shock 
velocity for the ordinary type Ib/c SNe than that deduced from the comparison between SN 
2003L and SN 2002ap and would, at the same time, further reduce the difference between 
"ordinary" and "engine- driven" type Ib/c SNe. 

The high shock velocities also suggest a high value for 7min- In an inverse Compton 
scenario, it may therefore be possible to observe the low energy c ut-off in the soft X-ray 



range. The rather hard X-ray spectrum (a ~ —1/2) observed by ISoderberg et al.l (120051 ) 
for SN 2003L could be caused by such an effect. Unfortunately, the low number of detected 
photons prevents a more detailed comparison. 

The observed radio spectra are steeper than expected from diffusive shock acceleration in 
a standard, non-relativistic shock. Inhomogeneities can influence several of the assumptions 
underlying the standard result. The fraction of electrons that scatters back across the shock 
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front as well as their average momentum gain depend on the scattering properties of the 
medium. The assumption of isotropic scattering may not be applicable for an inhomogeneous 
source structure. Furthermore, it is normally assumed that the velocity of the scattering 
centers downstream of the shock does not deviate too much from the average advection 
velocity. However, if the inhomogeneities are associated with regions of higher densities, 
they would be decelerated less at the shock front, i.e., have a higher velocity relative to the 
shock front than the rest of the medium. In cases when they take part in the scattering of 
the particles, the average momentum gain per shock crossing would go down, which would 
steepen the spectrum. 

In conclusion, the main points of the present paper can be summarized as follows: 

1) The presence of inhomogeneities in a synchrotron source will affect several of the 
main conclusions drawn from observations; for example, the observed correlation between 
total energy and shock velocity could, at least in part, be due to a varying covering factor. 

2) It is often hard to distinguish the effects of inhomogeneities from those of other 
mechanisms such as free-free absorption and cooling. It is argued that the broadening of the 
radio spectra and/or light curves is best understood as due to inhomogeneities; in particular, 
for most of the type Ib/c SNe in which cooling is likely to be negligible. 

3) Models, which attribute the X-ray emission in type Ib/c SNe to synchrotron radi- 
ation from a concave electron energy distribution produced in a non-linear shock, have a 
few implications that could limit their applicability; for example, the mass loss rate of the 
progenitor star needs to be substantially larger than deduced using a standard synchrotron 
model and the expected decline of the X-ray emission is slower than observed in many SNe. 

4) The observed correlation during the first few hundred days in type Ib/c SNe between 
the radio. X-ray and bolometric luminosities is the one expected from an inverse Compton 
scattering origin of the X-ray emission. An inhomogeneous source structure is consistent 
with equipartition conditions in the regions giving rise to the optically thin synchrotron 
emission. 

It is clear that inverse Compton scattering can only contribute significantly to the X-ray 
emission as long as the supernova is reasonably luminous. The X-ray emission during later 
phases needs to be attributed to another mechanism. Although free-free emission from a 
homogeneous source is not a viable option in man y cases, this would change with sufficiently 



large density inhomogeneities (e.g., IChugailll993l ). The radio emission is proportional to the 



product of the energy densities in the magnetic field and relativistic electrons. If both of 
these scale with the thermal energy density, one expects the radio emission to correlate with 
the free-free component of the X-ray emission. Hence, inhomogeneities due to a clumpy 
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medium would preferentially increase both of these emission components, while leaving the 
inverse Compton emission from either the thermal or relativistic electrons unchanged. 

The lack of broadening of radio spectra and/or light curves cannot be taken as evidence 
for a homogeneous, spherically symmetric source. For an unresolved source, assuming such 
a model results in a minimum value for the shock velocity. An independent indication for 
inhomogeneities is, therefore, a higher shock velocity deduced from observations at other 
wavelength ranges. When the shock velocity can be directly observed, unresolved inhomo- 
geneities will give rise to a radio brightness temperature lower th an expected. The low 



bright ness temperature in the type lib SN 1993 J was interpreted by iFransson &: Bjornsson 



(Il998l ) as due to cooling. Alternatively, it could be caused by inhomogeneities, which would 



then lower the deduced strength of the magnetic field. 



Thanks are due to an anonymous referee whose detailed and critical reading did much 
to improve the paper. 
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Fig. 1. — The observed values of the inverse Compton parameter {L^/ L\,o\Y / versus ra- 
dio luminosity L^.. The values have been normalized to those of SN 2003L. The lack of 
a correlation is evidence for an inverse Compton origin of the X-ray emission (see text). 
The values for 2003L, 19941 , 2001ig, 2002ap, and 2003bg are taken from th e compilation of 
Chevalier fc Fransson (|2006[). The other type Ib/c SN e shown are 2 008D flSoderberg et al.l 



2008h. 2004dk f|Wellons et al.l boid: IProut et al.lboilh . and 2009bb jSoderberg et al.lboiol : 
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Fig. 2. — The observed values of the inverse Compton parameter (Lx/Lboi)^/-^r versus time 
t after explosion for the same SNe as in Figure[Tl The weak variation with time constrains 
the properties of the inhomogeneities (see text). The downward pointing arrow for SN 19941 
corresponds to the upper limit of the X-ray emission on day 52. 



